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Abstract 

The steady state structure of an interface in an Ising system on a square lattice placed 
in a non-uniform external field, shows a commensurate -incommensurate transition driven 
by the velocity of the interface. The non-uniform field has a profile with a fixed shape 
which is designed to stabilize a fiat interface, and is translated with velocity v^. For small 
velocities the interface is stuck to the profile and is rippled with a periodicity which may 
be either commensurate or incommensurate with the lattice parameter of the square lattice. 
For a general orientation of the profile, the local slope of the interface locks in to one of 
infinitely many rational directions producing a devil's staircase structure. These "lock-in" or 
commensurate structures dissappear as Vg increases through a kinetics driven commensurate 
- incommensurate transition. For large fg the interface becomes detached from the field 
profile and coarsens with Kardar-Parisi-Zang exponents. The complete phase -diagram and 
the multifractal spectrum corresponding to these structures have been obtained numerically 
together with several analytic results concerning the dynamics of the rippled phases. Our 
work has technological implications in crystal growth and the production of surfaces with 
various desired surface morphologies. 
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I. INTRODUCTION 

Commensurate-in -commensurate (C-I) transitions^] have been extensively studied over 
almost half a century following early experiments on noble gases adsorbed on a crystalline 
substrate P] eg. Kr on graphite. Depending on coverage and temperature, adsorbates may 
show high density periodic structures the reciprocal lattice vectors (RLVs) of which are 
either a rational (commensurate) or irrational (in -commensurate) multiple of a substrate 
RLV. By changing external parameters (eg. temperature) one may induce phase-transitions 
between these structures. Recently, the upsurge of interest in the fabrication of nano- 
devices have meant a renewed interest in this field following a large number of experimental 
observations on "self-assembled" domain patterns (stripes or droplets) on epitaxially grown 
thin films for eg. Ag films on Ru(OOOl) or Cu-Pb films on Cu(lll)|3| etc. The alloy films 
often show composition modulations in the lateral direction forming patterned superlattices. 
These self-assembled surface patterns may have potential applications in the field of opto- 
electronics, hence the interest. In general, the whole area of surface structure modification 
has tremendous technological implications including, for example, the recording industry 
where magnetic properties are intimately connectedjj] to surface structure. 

Almost universally, C-I transitions may be understood using some version of the sim- 
ple Frenkel -Kontorova[5](FK) model, which models them as arising from a competition 
between the elastic energy associated with the distortion of the adsorbate lattice and sub- 
strate -adsorbate interactions. A complicated phase diagram involving an infinity of phases 
corresponding to various possible commensuration ratios (rational fractions) is obtained as a 
function of the two energy scales. In-between two commensurate structures one obtains re- 
gions where the periodicity of the adsorbate lattice is in -commensurate. All C-I transitions 
are equilibrium transitions in the sense that at any value of the relevant parameters, the 
structures observed optimize a free -energy. Indeed, despite its importance, the dynamical 
aspect of C-I transitions is a relatively unexplored domain. 

In this paper, on the other hand, we discuss a C-I transition entirely driven by kinetics. 
We show that (1) a simple Ising interface in a square lattice, held in place by a non-uniform 
external magnetic field, can have a variety of commensurate "phases" characterized by the 
local slope expressed in terms of the unit vectors of the underlying lattice and (2) it is possible 
to induce transitions in-between these phases by externally driving the interface with the 
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help of the field. The independent variables are therefore the average slope of the interface 
and its velocity both of which, as we show, can be externally controlled. Preliminary results 



from this work has been published elsewhere 



0,0. 



The dynamics of a 2-d Ising interface between the "up" and "down" spin phases at 
low teniperature T, in a (square) lattice driven by uniform external fields is a rather well 
studied^, 0] subject. The interface moves with a constant velocity, Voo, which depends on 
the applied field, h and the orientation 9 measured with respect to the underlying lattice. 
The interface is rough and coarsens with KPzQ] exponents a = 1/2 and (3=1/3 where a 
and f3 are the roughness and dynamical exponents respectively. We explore the possibility of 
driving such an interface with a pre-determined velocity Vf using an external non-uniform 
field which changes sign following a sharp sigmoidal profile forcibly stabilizing a stationary, 
macroscopically fiat interface at the region where the field crosses zero. We study systemat- 
ically the structure and dynamics of this "forced" Ising interface as the field profile is moved 
without change of shape at an externally controlable velocity Ve- We show that for low 
driving velocities Ve, the interface velocity Vf = Vg and the interface is stuck to the profile 
— the "stuck" phase. For larger Ve > Vod, the interface detaches. In the stuck phase the 
interface, though macroscopically fiat (i.e. a = /3 = 0) is patterned on the scale of the lattice 
spacing. It is these patterns which we show, undergo a series of C-I transitions determined 
by the Ve and the geometry characterized by the average slope of the interface in terms of 
the lattice vectors of the underlying square lattice. 

In the next section we introduce the model and briefiy sketch the main results from a 
mean field treatment. In section III we map our interface dynamics to the dynamics of a 
one dimensional "exclusion process" - a system of hard core particles on a line. In section 
III-A and B we present our results for the ground state structure and the dynamics within 
this model. Finally, in section IV we conclude. 

II. THE MODEL AND MEAN FIELD THEORY 

We consider here (see Fig.l) a one-dimensional interface Y{x,t) between phases with 
magnetization, 4>{x,y,t)> and (f){x,y,t) < 0, in a 2-d square lattice obeying single-spin 



fiip Glauber dynamics 



10| in the limit /i/ J, T/ J — > 0. Here J is the Ising exchange coupling 



and T the temperature. An external non-uniform field is applied such that h = hmax in 
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FIG. 1: (a) An Ising interface Y{x,t) (bold curved line) between regions of positive (marked +) 
and negative (marked — ) magnetization in an external, inhomogeneous field with a profile which 
is as shown (dashed line). The positions of the edge of the field profile and that of the interface 
are labelled Se and Sf respectively, (b) A portion of the interface in a square lattice showing 
a corner, (c) The interface velocity f/ as a function of the velocity of the dragging edge Ve for 
Ns = lOO(n), 1000(0), 10000(+) and p = 0.5. Ah the data (□,❖,+) collapse on the mean field 
solution (dashed line). Inset shows the graphical solution (circled) of the self-consistency equation 
for Vf] dashed line represents Vf = Vf. 

the +ve and —hmax in the -ve (f) regions separated by a sharp edge. The edge of the field 
(i.e. where the field changes sign) lies at S^- The front or interface, y(x, t), separates 
up and down spin phases. The interface is the bold curved line (Fig.l) with the average 
position Sf. When the edge is displaced with velocity Ve, the front, in response, travels with 
velocity vj. Parts of the front which leads (lags) the edge of the field experience a backward 
(forward) force pulling it towards the edge. The driving force therefore varies in both space 
and time and depends on the relative position of the front compared to that of the edge of 
the dragging field. In the low temperature limit the interface moves solely by random corner 
flips (Fig. ^b)), the fluctuations necessary for nucleating islands of the minority phase in 
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any region being absent. We study the behaviour of the front velocity vj and the structure 
of the interface as a function of Ve and orientation. 

Naively, one would expect fluctuations of the interfacial coordinate Y{x,t) to be com- 
pletely suppressed in the presence of a field profile. This expectation, as depicted by our 
main results (Fig^ FigI21and FigEj) is only partially true. While, as we show below, a mean 
field theory gives the exact behaviour of the front velocity f/ as a function of Ve (Fig^; 
small interfacial fiuctuations produce a dynamical phase diagram showing infinitely many 
dynamical phases (Fig El) and dynamic phase transitions (Fig. EJ. For Ve < foo the interface 
is stuck to the profile Vf = Ve- The stuck phase has a rich structure showing microscopic, 
"lock-in", commensurate ripples. These dissappear at high velocities through a dynamical 
commensurate- incommensurate (C-I) transition. 

Consider, first, a continuum coarse-grained description of the model shown in Fig. ^ . 
We assume that for /i/J, T/ J — » the magnetization is uniform everywhere except near the 
interface, Y{x,t) so that the magenetisation = 0(y — Y{x,t)). The field profile is given 
by h{y,t) = hmaxf{y,t) where f{y,t) = tanh((?/ — Vet)/x) and x is the width of the profile 
(see Fig^a)). Using Model A dynamics^ for and integrating out all degrees of freedom 
except those corresponding to the interfacial position, we obtain the effective dynamical 
equation satisfied by the interface. 



where Ai,A2 and A3 are parameters. Note that Eq. (j2.H) lacks Galilean invariance|ll| Y' 
Y + ex, x' ^ X — X2€t, t' ^ t. A mean field calculation amounts to taking Y = Y{t) i.e. 
neglecting spatial fiuctuations of the interface and noise. For large times (t —>■ 00), Y Vft, 
where f/ is obtained by solving the self-consistency equation; 



Vf = lim — Astanh ^ 



X 

= -A3sign(t;/ - Ve) (2.2) 

For small Ve the only solution to Eq. (j2.2|) is Vf = Ve and for Ve > Voo, where Voo = A3 
we get Vf = X3 = Voo- We thus have a sharp transition (Fig. E^c)) from a region where the 
interface is stuck to the edge to one where it moves with a constant velocity. How is this 
result altered by including spatial fiuctuations of y ? This question is best answered by 



mapping the interface problem to an assymmetric exclusion processjS], and studying the 
dynamics both analytically and numerically using computer simulations. 



III. FLUCTUATIONS: THE EXCLUSION PROCESS 

The mapping to the exclusion process follows by distributing Np particles among 



Any configuration of the system is specified by the set of integers {yi} where i/i denotes the 
location of the ith particle. In the interface picture i maps onto a horizontal coordinate {x in 
Fig-E)) and i/i as the local height Y{x). Each configuration {yi} defines a one-dimensional 
interface inclined to the horizontal with mean slope ta.n Of = l/p where p = Np/Ng. The 
yi satisfy the hard core constraint yj+i > ?/j + 1. The local slope near particle i is given 
by yi+i — yi and is equal to the inverse local density pi measured in a region around the z*^ 
particle. Alternatively, one associates a vertical bond with a particle and a horizontal bond 
with a holej^, in which case, again, we obtain an interface with a slope tan^^ = p/(l — p). 
The two mappings are distinct but equivalent. Periodic boundary conditions amount to 
setting yi+Np = 2/i ± Ng. Motion of the interface, by corner flips corresponds to the hopping 
of particles. In each time step {Np attempted hops with particles chosen randomly and 
sequentially^]), yi tends to increase (or decrease) by 1 with probability p (or g); it actually 
increses (or decreases) if and only if — yi > 1. The dynamics involving random sequential 
updates is known to indroduce the least amount of correlations among yi which enables 
one to derive exact analytic expressions for dynamical quantities using simple mean field 
argumentsjl^. The right and left jump probabilities p and q {p + q = I) themselves depend 
on the relative position of the interface yi and the edge of the field profile i/p + vj,. Note 
that this relative position is defined in a moving reference frame with velocity Vf{t), the 
instantaneous average particle velocity defined as the total number of particles moving right 
per time step. We use a bias Aj(t) = p — q = Asign(|/j — i/p — Vet) with A = 1 unless 
otherwise stated. In addition to the front velocity Vf, we also examine the behaviour of the 
average position, 



Ng sites of a 1-d lattice. The particles are labelled i = 1,2, 



Np sequentially at t = 0. 




(3.1) 



i=l,Np 
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FIG. 2: Variation of < y{t) > with t for Ve = 0.025 and p = 1.0. Lines denote analytic results 
while points denote Monte Carlo data for p = 1/5 (uppermost curve), 2/5 and an incommensurate 
p near 1/3. Inset (a)-(c) shows the corresponding ground state interfaces {yi — i/p). The arrows 
in (c) mark the positions of two discommensurations. 

and the width of the interface: 

^\t) = N;' Y1 < (y^(t)- < y.{t) >f > (3.2) 

i=l,Np 

as a function of time and system size A^^. Here, < yi{t) >= i/p + v^t. Angular brackets 



denotes an average over the realizations of 



hole symmetry for an exclusion processj^,^! is violated since exchanging particles and holes 
changes the relative position of the interface compared to the edge. 

We perform numerical simulations of the above model for Ns upto 10^ to obtain Vf for 
the steady state interface as a function of Ve as shown in Fig. ^^c). A sharp dynamical 
transition from an initially stuck interface with Vf = Ve to a free, detached interface with 
Vf = Voo = A(l — p) is clearly evident as predicted by mean field theory. The detached 
interface coarsens with KPZ exponentsy. Note that, even though the mean field solution 
for Vf{ve) neglects the fluctuations present in our simulation, it is exact. The detailed 
nature of the stuck phase {vf = and a bounded) is, on the other hand, considerably more 
complicated than the mean field assumption Y{x^t) = Y{t). 



:he random noise. Note that the usual particle 
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A. Ground state structure and the Devil's staircase 



The ground state of the interface in the presence of a stationary {ve = 0) field profile 
is obtained by minimizing E = 'Y^^{yi — i/ p — cY with respect to the set {y.i} and c. This 
maybe shown from Eq. (|2.H) by neglecting the terms non-linear in dYjdx] the resulting 
equation of motion, for small deviations of Y from the edge may be derived from the effective 
Hamiltonian E. The form of E leads to an infinite range, non-local, repulsive, interaction 
between particles in addition to hard core repulsion and the minimization is subject to the 
constraint that yi be an integer. For our system, the result for the energy may be obtained 
exactly for density p = m/n, an arbitrary rational fraction. For even m we have the following 
expression. 



E = - y [(A + c)2 + (A-c)2] + 

m ^-^ m m 

i=l,m/2 

i\ - cf + 

_ 1 1 1 1 1 1 

Where in the last equation we have minimized the expression with respect to c. 
Similarly for odd m we have the following expression for energy, 

E = - y [(A + c)2 + (A_c)2]+c2 

m ^-^ m m 

j=l,m/2-l 

= 1^(1 -A) (3-4) 

The resulting ground state profiles are shown in Fig. |21^ insets). The lower bound for 
E{p) is zero which is the energy for aX\ p = 1/n. For irrational p the energy is given by 
\im.m^^E{m/n) = 1/12 which constitutes an upper bound. For an arbitrary < p < 1 
the system {{yi}) therefore prefers to distort, conforming within local regions, to the nearest 
low-lying rational slope 1/p interspersed with "discommensurations" of density pa = \p — p\ 
and sign +ve (-ve) if these regions are shifted towards (away) from each other by 1. A plot 
of p(p) shows a "Devil's staircase" structure l|. In order to observe this in our simulation 
we analyze the instantaneous distribution of the local density of the particle -hole system to 
obtain weights for various simple rational fractions. A time average of these weights then 
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give us the most probable density p — distinct from the average p which is constrained to 
be the inverse slope of the profile. Increasing the width x of external field profile away 
from zero gradually washes out this Devil's staircase structure (Fig. O]). 
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FIG. 3: Devil's staircase structure in a plot of /5 Vs p for small velocities {v^ = 0.05) and 
X = 0.01, 1,5. 

As the velocity Ve is increased, steps corresponding to p = m/n (rational fractions) 
dissappear sequentially for Ve > 1/m so that for Ve > 1/2 only fractions of the form 1/n 
remain which persist upto Ve = foo- The locus of the discontinuities in the p{p) curve for 
various velocities Ve gives the phase diagram (Fig. |^ in the Vg — p plane. Note that the phase 
diagram for this C-I transition as given in an earlier publication 3| contained inaccuracies 
which have been now corrected. 
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FIG. 4: The dynamical phase diagram in Ve and p plane. The numbers on the p axis mark the 
fractions which determines the orientation of the lock-in phase. The three regions white, black 
and grey correspond to the rippled, the disordered and the detached phases respectively. 

B. Dynamics of the forced Interface 

For low velocities and density where correlation effects due to tlie liard core constraint are 
negligible, the dynamics of the interface may be obtained exactly. Under these circumstances 
the Np particle probability distribution for the ?/i's, P{yi,y2, ■ ■ ■ ,yNp) factorizes into single 
particle terms P{yi). Knowing the time development of P{yi) and the ground state structure 
the motion of the interface at subsequent times may be trivially computed as a sum of single 
particle motions. A single particle (with say index i) moves with the bias Aj(fet) which, in 
general, may change sign at y < i/ p + Vet < y + 1. Then P{yi) satisfies the following set of 
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master equations, 



P{yi) = -P{y,) + P{y, + 1) for > y + 1 

P{y,) = P{y, - 1) - P{y,) + P(i/, + 1) for y^ = y, y + 1 

P{y,) = -P{yi) + P{y, - 1) for yi < y. (3.5) 

The average position of the particle is given by < yiit) >= X]^=-oo -^(^/O 
spread by cr^(t) = X]^=-oo(2/*'~ < Vii^) Solving the appropriate set of master 

equations we obtain, for Ve « 1 the rather obvious steady state solution P{yi) = l/2{6y-^y + 
Sy-^y+i) and the particle oscillates between y and y + l. Subsequently, when i/ p + v^t > y + l, 
the particle jumps to the next position and P{yi) relaxes exponentially with a time constant 
r = 1 to it's new value with y ^ y + 1. For p = 1/n the entire interface moves as a 
single particle and the average position advances in steps with a periodicity of l/fg (see 
Fig. 121) In general, for rational p = m/n, the motion of the interface is composed of the 
independent motions of m particles each separated by a time lag of tl = 1/mVe- The result 
of the analytic calculation for small Vg and p has been compared to those from simulations 
in Fig. 121 for p = 1/5 and 2/5. For a general irrational p < 1/2, m -—>■ oo consequently, 

The forward motion of an irrational interface is accompanied by the motion of discom- 



mensurations along the interface with velocity Ve which constitutes a kinematic wave 
parallel to the interface. 

C. The C-I transition 

As the velocity fg is increased the system finds it increasingly difficult to maintain its 
ground state structure and for r > tl the instantaneous value of p begins to make excursions 
to other nearby low-lying fractions and eventually becomes free. Steps corresponding to 
p = m/n dissappear (i.e. p p) sequentially in order of decreasing m and the interface 
loses the ripples. The transition, as in the case of the FK model j^l] is characterized by well 
defined exponents. This may be seen by computing the spectrum of singularities 

/(a) = - 1)^,] - (g - 1)^, (3.6) 
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FIG. 5: The function p{p) for two different velocities fe/(l — p) = 0.1 and 0.5. Note that the steps 
corresponding to higher rational fractions tend to disappear with increasing Ve- Inset shows the 
corresponding multifractal spectrum /(a). Note that with increasing velocity, amin (see text). 

with a{q) = {d/dq){{q-l)Dg) and q = corresponding to the Devil's staircase in 

p(p). Dq is the Hausdorff dimension and Dg the generalized dimensions. Dg's are obtained 
by solving for 

Y.ipl/lf^''^) = 1. (3.7) 

i 

The changes in p determined the scales /j of the partition (defined following the Farey 
construction) whereas the changes in p were defined to be the measu res k . The high -order 
gaps in the vicinity of primary steps (the 1/n fractions) scale like|l3| pi ~ z"™*" where 
o^min = -Doo- Near these steps p ^ {p — Pmax)^ where Pmax is the maximum value of p for a 



step. This universal critical exponent |l4j| ^ = amin has been determined to be 0.71 ± .001 
from our data at fe/(l — p) = 0.1. The exponent ^ determines the stability of the rippled 
pattern to small changes in the orientation of the external field (p). As increases, ^ — > 

(Fig. ED. 

It is important to realize that the C-I transition seen in our system is driven by fluc- 
tuations of the local slope and therefore different from the C-I transition in a mechanistic 
Frenkel -Kontorova[l] model appropriate for domain structures arising from atomic misfits. 
The non-local energy E{p) and the non-linear constraint Ui = integer makes it extremely 
difficult to devise a natural mapping of this problem onto an effective F-K model. However, 
an approach based on the Langevin dynamics of particle of a single particle with coordinate 
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p' diffusing on a energy surface ^ven by, F{p') = E{p') + K,{p' — pY can obtain the main 
qualitative results of this section j7]. Here, k is an arbitrary constant which ensures that the 
long time limit of p' = p. 



IV. CONCLUSION 



In this paper, we have studied the static and dynamic properties of an Ising interface in 
2-d subject to a non-uniform, time- dependent external magnetic field. The system has a rich 
dynamical structure with infinitely many steady states. The nature of these steady states 
and their detailed dynamics depend on the orientation of the interface and the velocity of 
the external field profile. In future we would like to study the statics and dynamics of such 
forced interfaces in more realistic systems eg. a liquid-solid interface produced by coupling 
to a patterned substrate |3] . The authors wish to thank S. S. Manna for discussions. A.C 
acknowledges C.S.I.R. Govt, of India for a fellowship and S.S. thanks the Department of 
Science and Technology, Govt, of India, for financial support. 
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